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In the global fit of the SM using Gfitter, electroweak precision observables as well as constraints 
from direct Higgs searches have been compared with state-of-the-art electroweak predictions. 
We use the most recent results for direct Higgs searches at LEP and Tevatron and the latest 
measurements of mt and Mw- Example results are an estimation of the mass of the Higgs 
boson (Mh = 116.3 li 5 3 6 GeV) and a forth-order result for the strong coupling constant 
(as (Ml) = 0.1193 ± 0.0028(exp) ± O.OOOl(theo)). A fit of the oblique parameters (STU) to 
the electroweak data is performed, in order to constrain physics beyond the Standard Model. 
For instance, the parameter space of the Littlest Higgs Model with T-parity can be restricted 
via the oblique parameters. In addition, fit results for a model with an extended Higgs sector 
(2HDM) using mainly observables from the B and K physics are presented. 

1 Introduction 

Precision measurements allow us to probe physics at much higher energy scales than the masses 
of the particles directly involved in experimental reactions by exploiting contributions from 
quantum loops. Prominent examples are the electroweak precision measurements, which are 
used in conjunction with the Standard Model (SM) to predict via multidimensional parameter 
fits unmeasured quantities like the Higgs mass. 

Several theoretical libraries within and beyond the SM have been developed in the past 
containing the pertubative calculations of the SM and new physics models for the electroweak 
observables. However, most of these programs are relatively old, were implemented in outdated 
programming languages, and are difficult to maintain with respect to the theoretical and exper- 
imental progress expected during the forthcoming era of the LHC. These considerations led to 
development of the generic fitting package Gfitter [1] , designed to provide a modular framework 
for complex fitting tasks in high-energy physics. Gfitter is implemented in C++ and relies on 
ROOT functionality. The package allows a consistent treatment of statistical, systematic and 
theoretical errors, possible correlations and inter-parameter dependencies. 
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In this paper we present Gfitter results of the global SM fit to the electroweak observables 
as well as an estimate of the oblique parameters, which can be also used to constrain the 
parameter space of the Littlest Higgs Model (LHM). In addition, a fit of a Two Higgs Doublets 
Model (2HDM) is performed using B and K physics observables. 

2 The Global SM Fit 

In the global electroweak fit the state-of-the-art calculations of the electroweak precision observ- 
ables are compared with the most recent experimental data to constrain the free parameters of 
the fit and to test the goodness-of-fit. The free parameters of the SM relevant for the global elec- 
troweak analysis are the coupling constant of the electromagnetic, weak, and strong interactions, 
as well as the masses of the elementary fermions and bosons. Due to electroweak unification and 
simplifications arising from fixing parameters with insignificant uncertainties compared to the 
sensitivity of the fit, the number of free fit parameters can be reduced. The remaining floating 
parameters in the fit are the coupling parameters Aa[^ d (M|) and «s(M|), the masses Mz, rn c , 
rn b , m top , and M H . 

In Gfitter a complete new library of the electroweak precision observables as measured by 
the LEP, SLC, and Tevatron experiments has been implemented. State-of-the-art predictions in 
the one-mass-shell scheme are used. In particular, the full two-loop and leading beyond two-loop 
corrections are available for the predictions of M\y [2] and sin 2 6 l eS |3|4j . The implementation 
of the NNNLO pertubative calculation of the massless QCD Adler function [5j, contributing 
to the vector and axial radiator functions in the prediction of the Z hadronic width, allows 
to fit the strong coupling constant with a unique theoretical accuracy. Wherever possible the 
calculations have been cross-checked against the ZFITTER package [6]. More details on the 
theoretical computations in Gfitter can be found in pQ. 

The following experimental measurements are used: The mass and width of the Z boson, the 
hadronic pole cross section o~^ ad , the partial widths ratio R®, and the forward-backward asym- 
metries for leptons A^, have been determined by fits to the Z line-shape measured precisely 
at LEP (see [7J and references therein). Measurements of the r polarization at LEP [7J and the 
left-right asymmetry at SLC [7J have been used to determine the lepton asymmetry parameter 
Ai . The corresponding c and 6-quark asymmetry parameters A c ^ , the forward-backward asym- 
metries A™ , and the widths ratios R® and R®, have been measured at LEP and SLC [7j. In 
addition, the forward-backward charge asymmetry measurement in inclusive hadronic events at 
LEP was used to directly determine sin 2 9 l s |7J . For the running quark masses m c and mj the 

world average values are used. For Aa^ d (Af|) we take the phenomenological result [8]. For the 
W width we use the official combined LEP and Tevatron result, while for the W mass we take 
also into account the recent DO measurement [9] leading to our private preliminary combined 
value of M\y = (80.399 ± 0.023) GeV. In case of the top mass the latest combined result [10] 
mt = (173.1 ± 1.3) GeV is used. The direct searches for the SM Higgs Boson at LEP [11] and 
the most recent results from the Tevatron [12) . leading to a 95% confidence level (CL) exclusion 
for Mh < 114.4 GeV and at Mh = [160, 170] GeV respectively, are included using a Gaussian 
approach that quantifies the difference between the observed test statistics (the log-likelihood 
ratios) and the expected values for the s+b hypothesis using the values of the respective con- 
fidence level (CLs+b)- A contribution to the x 2 estimator of the fit is derived for each Higgs 
mass. More details of the procedure can be found in pQ. We perform global fits in two versions: 
the standard ("blue-band") fit makes use of all the available information except for the direct 
Higgs searches and the complete fit uses also the constraints from the direct Higgs searches. 

The standard (complete) fit converges at the global minimum value Xmm = 16.4 (Xmm = 
17.9) for 13 (14) degrees of freedom, corresponding to a p-value of 0.228 ± 0.004_o.oo2 (0.204 ± 
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Fig. 1: Ax 2 as a function of Mh for the complete fit. 
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Fig. 2: Determination of Mh excluding all the sensi- 
tive observables from the standard fit, except for the 
one given. 



0. 004_o.oo2) derived from MC toy experiments. The estimation for Mh from the standard fit, 

1. e., without the direct Higgs searches is Mjj = 83I23 GeV and the 2a and 3a intervals are 
respectively [42, 158] GeV and [28, 211] GeV. The complete fit represents the most accurate 
estimation of Mh considering all available data. The resulting Ax 2 curve versus Mh is shown 
in Figure [TJ The shaded band indicates the influence of theoretical uncertainties, which are 
included in the fit with a flat likelihood within the allowed ranges. The inclusion of the direct 
Higgs search results from LEP leads to a strong rise of the Ax 2 curve below Mjj = 115 GeV. The 
data points from the searches at the Tevatron, available in the range 110 GeV < Mh < 200 GeV 
increases the Ax 2 estimator for Higgs masses above 150 GeV beyond that obtained from the 
standard fit. The estimation for Mh from the complete fit results in Mh = 116.3 Ij 5 ^ 6 GeV 
and the la interval is reduced to [114, 145] GeV. 

In figure [2] only the observable indicated in a given row of the plot is included in the fit. Only 
for the four observables providing the strongest constraint on Mh, namely Ai(LEP), ^(SLD), 
j4p B and My/, the Higgs mass is determined. The compatibility among these measurements can 
be estimated by repeating the global fit where the least compatible of the measurements (here 
Ap B ) is removed, and by comparing the Xmin estimator obtained in that fit to the one of the full 
fit (here the standard fit). To assign a probability to the observation, the Ax^ in obtained this 
way must be gauged with toy MC experiments to take into account the "look-elsewhere" effect 
introduced by the explicit selection of the pull outlier. We find that in (1.4 ± 0.1)% ("2.5<r") of 
the toy experiments, the Ax^m found exceeds the Ax^ in = 8.0 observed in current data. 

The strong coupling at the Z-mass scale is determined by the complete fit to as (M%) = 
0.1193±0.0028±0.0001 where the first error is experimental and the second due to the truncation 
of the pertubative QCD series. 

Figure [3] compares the direct measurements of M\y and nit, shown by the shaded/green 
la bands, with the 68%, 95%, and 99% CL obtained for three sets of fits. The largest/blue 
(narrower /yellow) allowed regions are derived from the standard fit (complete fit) excluding the 
measured values in the fits. The inclusion of the LEP and Tevatron Higgs searches significantly 
impacts the constraints obtained. Figure [3] allows to compare the indirect and direct determina- 
tion of the Mw and nit. So far the indirect determinations and the direct measurements are in 
good agreement. The third set of fits (narrowest/green) results from the complete fit including 
the measured values. Hence, it uses all available information and leads to the narrowest allowed 
region. 
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Fig. 3: Contours of 68%, 95%, and 99% CL obtained 
from scans of fits with fixed variable pairs Mw vs. m< 
for three sets of fits. The horizontal bands indicate 
the la regions of measurements (world averages). 
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Fig. 4: Allowed contours of 68% and 95% CL in the 
(T,S)-plane obtained from fits with m t = 172.4 GeV 
and M H =116, 350, and 1000 GeV. 



3 Oblique Parameters 

A common approach to constrain physics beyond the SM using the global electroweak fit is 
the introduction of oblique parameters, which assumes that the contributions of new physics 
models only appear through vacuum polarization. Most of the effects on electroweak precision 
observables can be parametrized by three gauge self-energy parameters (S, T, U) introduced by 
Peskin and Takeuchi [13J. S (S + U) describes new physics contributions to neutral (charged) 
current processes at different energy scales, while T measures the difference between the new 
physics contributions of neutral and charged current processes at low energies {i.e., T is sensitive 
to isospin violation). Further generalizations like additional corrections to Zbb couplings [2] 
can be also taken into account. 

The constraints on the STU parameters are derived from the fit to the electroweak precision 
data, presented in section [2j The STU parameters replace Mh and mt as free parameters in 
the fit. The following fit results are determined from a fit assuming as(M|) = 0.1193 db 0.0028, 
m t = 172.4 GeV, and M H = 116 GeV (in parentheses M H = 350 GeV): 

S = 0.02 (-0.06) ± 0.11 

T = 0.05 ( 0.15) ± 0.12 (1) 
U = 0.07 ( 0.08) ± 0.12 

Since U is generally small in new physics models and only constraints by the mass and width of 
the W boson, the STU parameter space is often projected to a two-dimensional parameter space 
in which the experimental constraints are easy to visualize. Figure [4] shows the 68%, 95%, and 
99% CL allowed contours in the (T,5)-plane for three different assumptions for Mh- In any case 
the oblique parameters are small, i.e., possible new physics models may affect the electroweak 
observables only weakly. 



4 Littlest Higgs Model with T-Parity 

The fine-tuning problem of the Higgs mass parameter (hierarchy problem) is one of the driving 
arguments to consider physics beyond the SM. Besides supersymmetric extensions of the SM, 
Little Higgs theories provide a way to tackle the hierarchy problem. The generic structure of 
Little Higgs Models is a global symmetry broken at a scale / (around 1 TeV) where new gauge 
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Fig. 5: Allowed contours of 68%, 95%, and 99% CL 
obtained from scans of fits with fixed variable pairs / 
and Mb {s\ = 0.45). The parameter F is a quantita- 
tive measure of fine-tuning (larger values of F corre- 
spond to larger degree of fine-tuning). 




Fig. 6: Simple overlay of the 95% CL exclusion regions 
in the (tan /3, M H ± )-plane from individual 2HDM con- 
straints and the toy-MC-based result (black solid line) 
from the combined fit overlaid. 



bosons, scalars, and fermions exist canceling the one-loop quadratic divergences to the Higgs 
mass from the SM particles. 

The Littlest Higgs Model (LHM) [15] is based on a non-linear la model describing an 
SU(5)/SO(5) symmetry breaking. T-parity conservation can provide a possible candidate for a 
dark matter WIMP (similar to R-parity conservation in supersymmetry). In addition T-parity 
forbids tree-level contribution from heavy gauge bosons to the electroweak precision observables. 
The STU parameters of the oblique parameter fit are replaced by the calculations of the cor- 
responding ones in the LHM [16J. The new floating parameters of the fit are: / the symmetry 
breaking scale, s\ ~ mT-/mT+ in leading order the ratio of masses of the T-odd and the T-even 
state from the LHM top sector, and 5 C a order one-coefficient, which exact value depends on 
detail of UV physics. The latter parameter is treated as theory uncertainty in the fit S c = —5. ..5. 

Figure [5] shows the 68%, 95%, and 99% CL allowed contours in the (M#,/)-plane for a fixed 
value of s\ = 0.45. Contributions of the T-odd partners of light fermions to the T parameter are 
neglected. This assumption is justified as long as the T-odd fermions are sufficiently light. The 
parameter F is a quantitative measure of fine-tuning, indicated by the black lines. Since larger 
values of F correspond to larger degree of fine-tuning, large values of Mh are more preferred than 
small values. Therefore, large Higgs masses are not only allowed by the electroweak precision 
data, but they are also favored in terms of fine-tuning. 

5 Two Higgs Doublet Model 

In the Type-II 2HDM, we constrain the mass of the charged Higgs and the ratio of the vacuum 
expectation values of the two Higgs doublets using current measurements of observables from the 
B and K physics sectors and their most recent theoretical 2HDM predictions, namely R® |7|17| . 
the branching ratio (BR) of B — ► X s ^y |18|19j . the BR of leptonic decays of charged pseudoscalar 
mesons (B — ► tv |20|21j . B — > fiv |22|21j and K — > jsu [23J) and the BR of the semileptonic 
decay B -> Dtv |24l25] . 

For each observable, individual constraints have been derived in the (tan/3, Mg±) plane. 
Figure [6] displays the resulting 95% excluded regions derived assuming Gaussian behavior of the 
test statistics, and one degree of freedom. The figure shows that R? is mainly sensitive to tan j3 
excluding small values. BR(1? — > X s j) is only sensitive to tan/3 for values below ~1. For larger 
values it provides an almost constant exclusion of a charged Higgs lighter than ~260 GeV. For 
all leptonic observables the 2HDM contribution can be either positive or negative since signed 



terms enter the prediction of the BRs resulting in a two-fold ambiguity in the (tan/3, m#±) 
space. 

In addition, we have performed a global fit combining the information from all observables. 
For the CL calculation in the two-dimensional plane we performed toy MC tests in each scan 
point which allows to avoid the problem of ambiguities in the effective number of degrees of 
freedom. The 95% CL excluded region obtained are indicated in Figure [6] by the area below 
the single solid black line. We can exclude a charged Higgs mass below 240 GeV independently 
of tan/3. This limit increases towards larger tan/3, e.g., M H ± < 780 GeV are excluded for 
tan/3 = 70. 
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